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In the present paper, we analyze the consequences of scaling hypotheses on dynamic functions, 
as two times correlations C(t,t'). We show, under general conditions, that C(t,t') must obey the 
following scaling behavior C(t,t') = (f>\{t)^^S{f3), where the scaling variable is = /3(<f>i(t') / fa(t)) 
and 4>i(t'), 4>i{t) two undetermined functions. The presence of a non constant exponent f((3) signals 
the appearance of multiscaling properties in the dynamics. 



I. INTRODUCTION 



The introduction of scaling concepts to describe equilibrium and off equilibrium dynamics in statistical mechanics 
was originally motivated by experimental and simulation data about, for instance, structure factor, pair correlation 
functions, response functions. Actually, the study of several classes of materials with complex dynamical properties 
as magnets, polymers, glasses, and several other thermal systems, and even non-thermal systems as granular media, 
, has shown the presence of some general scaling features In order to formulate a coherent scaling approach 

Q\ • to the dynamics of systems out of equilibrium, in this paper we resort to a general scheme developed in 1971 [0. 
This approach also reproduces as a particular case the multifractal and the multiscaling formalism, which has been 
applied to a large variety of phenomena such as turbulence, random resistor networks, self organized criticality, 
spinodal decomposition and many more [j5|-]l3f . The general scaling formulation applied to systems out of equilibrium 
stems from the hypothesis of invariance of two time functions, as autocorrelation functions, under a general scaling 
' transformation with the only requirement that the transformation obeys group properties. 

For defmiteness let's consider a two time correlation function, C(t,t'), which for example could be the density- 
density autocorrelation function in a supercooled liquid or the spin-spin correlation function in a magnetic system. 
We suppose that the system is prepared at time t = and is probed at two subsequent time t' and t. When the 
I ■ system is out of equilibrium, C(t, t') generally depends on both t' and t. Whenever the relaxation characteristic time, 
r, is very large or infinite, we may make the following asymptotic scaling ansatz valid for t and t' large but smaller 
than t: by rescaling the system lengths by a factor I, if t and t' are opportunely rescaled, we may expect that the 
autocorrelation function scales as £-f(*>* ), where the exponent f(t,t') is in general dependent on t and t' To be 
precise we assume that the function C(t,t') has the following general scaling property 



c(t,t f ) = r^ t ' t ">c(t,t') (i) 



' under a general time rescaling "non mixing" transformation, which satisfies group rules, as 

t = F 1 (t,l) t' = F 2 (t',l) (2) 

with the condition that Fi(x, 1) = x (i = 1,2). The above transformations are "non mixing" in the sense that F\ 
(resp. F2) depends only on t (resp. t'). The requirement that the transformation obeys group properties imposes 
some constraints on the functions C(t,t') and f(t,t'). Interestingly, under these assumptions, we find that C(t,t') 
can be synthetically expressed in the following way: 

C(t,1?) = Mt) m S(fi), (3) 
where the scaling variable, (3, has the following form: 

[3 = fait')/ fait) . (4) 

Here the fa (i = 1,2) are two unknown functions fixed by the transformations given in eq. (^]). Eq. (Q) in the 
particular case f(f3) — was obtained in Ref. |lj| using different arguments. Notice that, whenever / is not a 
constant, a "multiscaling" dynamical behavior is found in the dynamics, an interesting issue to check in models as 
well as experiments and simulations. 
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II. THE GENERAL "NON MIXING" CASE 

In what follows we give a demonstration of what summarized above. As shown in Ref. 0], the general transforma- 
tions of eq. (|^) implies that exist a couple of functions, fa(t) and 4> 2 (t'), which under rescaling exhibits the following 
property: 

Mt) = Mt)/i fa(t') = Mt')/i (5) 

These equations state that the "true" scaling variables are the (pi's, and that, whenever the functions fa's are 
invertible, eq.s (||) can be expressed in the following way: 

t = F 1 (t,l) = fa^ 1 (fa(t)/l) i> = F 2 (t',l) = fa^ 1 (fa(t')/l) (6) 

where fa~ x is the inverse function of fa (i.e., 4'Y 1 (fa(x)) = x). 

Let's now study how the above group properties influence the structure of the function C(t,t'). The group rules 
impose that if we scale t and t' by a factor l\ and later by a factor l 2 this should be equivalent to rescale them by a 
factor li ■ l 2 . More formally, we can express this condition as: 

C(F 1 (F 1 (t ) l 1 ) ) l 2 ),F 2 (F 2 (t',h),l 2 )) = (hl 2 )-^C(t,l/) . (7) 
Substituting eq.s (|j) and (|2|) in the above relation, one is led to a simple equation which states that: 

±f(t(l),t'(l))^0 (8) 

where, by definition, t(l) = F\{t,l) and t'(l) — F 2 (t,l). By inserting eq.(^) in eq. (||), one finds that 
/(0rWO^ 2 ~WO) = fifa: 1 (fa/1), fa^ 1 (fa/1)), i-e., we have that f(faZ l {fa/l)A 2 \fa/l)) = /(*,*')• Now, 
by taking I = fa we obtain that /(*,*') = f Ui 1 (1) , fa^ 1 ((f> 2 / fa)) , that is to say f(t,if) = f(fa(t')/fa(t)). 

Analogously, by inserting eq.(||) in eq. (|l|), and choosing I — fa, we find the scaling form for C(t,t') that we 
anticipated in eq. (|3j) above. In such a way we also individuate the scaling function S: Six) = C (fa^ 1 (1) , c/) 2 1 (x)). 

Thus we proved that in presence of scaling properties as those written in eq.s (Q) and (||), the asymptotic functional 
form of the scaling of C(t,t') is characterized by the asymptotic behavior of the "true" scaling variables </>i and fa, 
as written in the general result of eq. (^) . 

For sake of clarity we have dealt with a two variables function, C(t,t'), but analogous properties may be proven 
for a many variables function, C(t\, t 2 , t n ). In this case, if the generic variable undergoes a scale transformation 
U = Fi(U,l) (i e {1, ...,n}), we have: 

C{h,t 2 ,...,t n ) = fa(t i y^>-^S(p 2 ,. ..,/?„) (9) 

with (i > 1) ^ = fa(ti)/fa(ti) As before, /(xi, ...,x n -\), and S{x\, ...,x n -\) are undetermined functions. 



A. Some examples 

In many physical cases we might generally expect that the two times, t and t' , scale in the same way, i.e., fa — fa = fa 
Below we explicitly list few interesting examples in this category. 

A simple situation corresponds to a scaling function, fat), which is asymptotically a power law in t (see ref.s in dJ), 
i.e., one has fat) ~ t x l z , and the scaling of C(t,t') is: C(t,t') = t^ 1 '^S(t'/t). In most cases we expect the exponent 
/ to be a constant, so that: 

C(t,t') =t f ' z S{t'/t) . (10) 



Asymptotically, the "power scaling" of eq. (|1C|) , is found in several toy models for glasses as in a "phase space" model 
Jl6| , the Backgammon entropic barriers model jH| , the Queens long range interactions model |jl9| , in solvable models 
of interacting particles in high dimensionality [2^, or in a kinetically constrained lattice gas |2l| (see also references 
in ||). Several of these cases are characterized by / — 0. But in general one might expect cases with / different 
from 0. For instance in the Bak-Sneppcn SOC model, the two time function P(t,t'), describing the return of activity 
to a site at time t which was most recently active at time t' , for an avalanche started at t — 0, seems to have a 
scaling of the form: P(t, t') = t ZBS V(t! jt) p2fl, with constant zbs- In a model of direct polymers in a random media 
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a similar behavior is found for the off equilibrium "overlap function": q(t,t') = t~ x q(t'/t) (23). In some models of 
non-linear diffusion equations correlation functions have been shown to have a "power law" scaling structure 
of eq. ( |Io| ) with constant exponents / and with a scaling function S which is itself a power law. In the framework 
of out equilibrium dynamics, phenomena as coarsening or, more generally, phase ordering kinetics in "standard" 
Ginzburg-Landau magnets usually show correlation functions which are asymptotically characterized (see [||) by the 
above scaling of eq. (|f 0|), which is often called "simple or full or naive aging". In some discussions of glassy relaxation 
also a more complex, "interrupted aging", scenario was proposed pUTq], in which the long time regime of the two time 
autocorrelation function scales as C(t,t') — S(t' /t 1+fJ -). In the present picture this corresponds to two different power 
exponents for faft) ~ t x / z and <j) 2 {t') ~ t'V*+M. 

In the case where / = f(t'/t) is a non-trivial function of the scaling variable /3 = t'/t, one finds a multiscaling 
dynamical behavior. This is analogous to the multiscaling found, in different context, by Coniglio and Zannetti in the 
spinodal decomposition of the N = oo Ginzburg Landau model with conserved order parameter or the one proposed 
also for the density profile of the DLA model [|l3| . 

In the previous cases, the scaling variable was the ratio of powers of the two involved times, however in different 
situations, as for instance in the limit in which the exponent 1/z goes to zero, one may expect to have a logarithmic 
behavior for <j>: <j>{t) ~ ln(t). This situation gives as scaling structure: C(t,t') = ln(t) / ( ln ( t ')/ ln ( t »5(ln(t')/ ln(i))- In 
many cases one has / = 0, namely: 

C(M') = S(m(t')/ln(*)) . (If) 

The "logarithmic scaling" of eq. ( |TT| ) is found in several systems. An example of diffusion which shows the "logarithmic 
scaling" is the one dimensional Sinai's model with a random local bias. In this case, for instance, the two time residency 
probability asymptotically has a scaling form given by eq.(|ff|) with a scaling function S(/3) which is an exponential 
corrected by a power law in f3 = ln(t')/ln(£) p5fl . Random Field systems also show "logarithmic scaling" [||,[|, 
but also experimental random exchange Ising ferromagnets [^6| , among many other belong to this category. 

Logarithmic kinetics have been also recently experimentally observed in the amorphous-amorphous transformations 
in some glasses under high pressure |27j| . Interestingly, also non-thermal systems as granular media, shaken at low 
vibration amplitudes, present a non-trivial out of equilibrium dynamics, where numerical calculations on different 
models jl7j suggest a "logarithmic" scaling in the relaxation of the two time density correlation function as in eq.(^fj). 

The scenario about other disordered systems such as spin glass models is still controversial. To describe numerical 
calculations and to fit experimental data of relaxation in the thermoremanent magnetization of some spin glasses 
several proposals as power scaling eq.(|IO|) or logarithmic scaling eq.(|ff]) have been made p^-|3f|^[|. Also in recent 
computer simulations of a Lennard- Jones off equilibrium glass model the asymptotic behavior of the autocorrelation 
function was suggested to have a "logarithmic" scaling ||32|| as opposite to a "power law" scaling previously proposed 

III. THE "MIXING" CASE 

Up to now we have dealt with "non mixing" scale transformations, as in eq. where the scaling of each of the 
variables doesn't depend on the other. However, situations where "mixing" is present might be possible. Formally the 
case of mixing may be dealt with as the non-mixing one, however the results are too general to be of immediate practical 
use. For sake of completeness, we just show them. In the mixing case one finds that eq. (^) must be substituted 
by: C(t,t') — 4>i{t,t')^^S{(3), where the scaling variable, /?, is now: (3 — fait, t')/4>i{t, t'). Here, as before, the fa 
(i = 1,2) are two unknown functions fixed by the mixing transformations t — Fi(t,t',l) and V = F2(t,t',l) (with 
Fi(t,t', 1) = t and F2(t,t', 1) = t 1 ). The above result may be of scarce use because any function of two variables, 
C(t,t'), may be written as above in terms of other two functions, fa{t,t') and fa(t,t'). 

However, it may be interesting to work out a specific example of mixing transformations, which shows how one may 
recover, from simple scale principles, a multifractal scaling structure. 

A. An example of mixing 

While the function C(t,t') has the general scaling property of eq. ([!]), we now suppose that the rescaling transfor- 
mations of t and t' have the following specific form under a scale change of extension I: 

t = t/l t' = t'/l z ^ . (12) 



3 



Interestingly, within this context, we find that the scaling of C(t,t') is restricted to have the following structure: 

C(t, t') = t f{p) S{(3) . (13) 
Here the scaling variable, [3, has only two possible forms: either it is a ratio of powers of the two times 

P = t'/t z (14) 
with z = const (corresponding to a non- mixing case previously described), or 

where H((3) is an undetermined function. 

The scaling form (|T4|) corresponds to the case z(t, t') — const, which is one of the non-mixing cases we dealt with 



before. The scaling form given in eq. (15) corresponds instead to a non constant scaling exponent z — z(t,t') in 
which thus gives a mixing transformation of t and t' . Actually, it turns out that the only possible solution for a non 
constant z is z(t,t') = with [3 given in eq. (fi"5"|). In this case the scaling variable is asymptotically logarithmic in 
the two times, (3 = ln(t')/ln(t) + 0(l/ln(t)). This kind of scaling for different variables was proposed, for instance, 
for the multifractal description in turbulence in the DLA model [Q, in Self-Organized-Critical (SOC) models 

H or in voltage distribution of random resistor networks These scaling forms, differently from ordinary critical 

phenomena, are characterized by a continuity of scaling exponents. 

For definiteness it is interesting to work out the simple case where the function H{(3) is linear in (3, situation which 
might generically correspond to the case of very long times t and small (3. By writing H({3) = (3 ln(fo) — hi(ig) (where 



tg and t' are constants), from eq. (15) one finds that (3 — ln(t'/t' )/ ln(t/io). This case corresponds for instance to the 
multifractal scaling proposed by Kadanoff et al. in Ref. || to describe the avalanches size distribution in the context 
of SOC models. 

Below we work out the example of mixing transformation of eq. (|l2j ) in details. As in the general case above, we 
have to impose the group rules on the scale transformations. For the transformation of the variable t and t' this 
implies (i = 1,2): 

F i (F 1 (t,t',l 1 ),F 2 (t,t',l 1 ),l 2 )=F i {t,t',h-l2) . (16) 
For a transformation as in eq. (^|), this assertion simply imposes that: 

j lZ (t(l),t'(l)) = (17) 



The above eq. (17) has two kind of solution. The first is the trivial one: z = const.. The second is non-trivial and 



has the following form: 

_ ln(t') H(z) 
z ~ Htj + Htj 



(18) 



where H(z) is a generic function. The latter may be obtained as follows. Since the function z(t,t') = z(t/l,t' /l z ) 
is invariant under rescaling eq. (|l7|), we can write that z(t/l,t' /l z ) — z(t,t'). By fixing I = t, we obtain : z{t,t') = 
z(l,t'/t z ). Thus z is a function of the single variable t' /t z , and we can write z(t,t') = g(t' /t z ). Here we have defined 

g(t'/t*)=z(l,t'/t z ) . (19) 

By inverting the above relation we have: t'/t z — g^ 1 (z), and passing to the logarithms we recover z = ln(i')/ln(£) + 
H(z)/ln(t), where we have introduced the unknown generic function H(z) = — ln(g _1 (z)). Thus we have found the 



solution given in eq. (18). This result states that, for fixed z, whenever ln(t') and ln(£) are large enough, we have 
z = \n(t')/Ht)+0(lJ\n(t)). 

By imposing then group properties to the function C(t, t') itself (see eq. (0)), we obtain eq. (|8|), which, after insertion 
of eq. ©, implies that f(t/l,t'/l z ) = f(t,t'), and by taking I = t, as above, we have that f(t,f) = /(l, t'/t^^). 

Whenever z is a constant we recover a non- mixing case described in the previous section. Let's now suppose that 
z is not a constant and is given by eq. (|l8|). In this case, we can prove that the exponent f(t,t') is a function of 
the single variable z: f(t,t') — f(z). In fact, as before we have that f(t,t') — f(l,t'/t z ) = f(l,g^ 1 (z)), i.e., / is 
a function of the variable z. Analogously one proves that C(t,t') = t^ z >S(z), where the scaling function S is now 
S(x) = C(l, g~ 1 (x)). From this result, in the asymptotic limit of large t' and t, we recover the scaling form for C(t, t'), 
given m eq.s (||) and ©. 
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IV. CONCLUSIONS 



We expect that the present approach may be useful to describe general properties of dynamical functions in physical 
systems when their characteristic times diverge, since in such a situation, like close to usual critical points, scale 
invariance should be reasonably present. Actually, one observes diverging characteristic times typically when out 
of equilibrium dynamics phenomena become important, i.e., when an explicit dependence of functions as C(t,t') on 
both times (and not on their difference) is observed, a fact which is some time generically called "aging" . In this 
perspective the structural properties of scaling described here should be naturally associated to out of equilibrium 
dynamics (i.e., to "aging") effects. Interestingly, we have pointed out that in a broad variety of physical systems, 
ranging from magnets, to polymers, glasses, or spin glasses, random fields, random ferromagnets, granular materials, 
diffusive systems, etc., one observes scaling properties of dynamical functions which may be well inserted in the 
framework reported above. 

We have shown in full generality that a generalized homogeneous function, C(t,t'), which acts as in eq. ([j]) under 
the scale transformation of its variables given in eq. (g), must obey the scaling behavior of eq. (||). In this theoretical 
framework, also a multiscaling or multifractal behavior is admissible in the dynamics. It would be interesting to 
understand if it exists in real dynamical systems. 

The present approach is not restricted to scaling of dynamical functions. We have seen that it is relevant to 
describe, as well known, usual scaling in standard critical phenomena, but it also describes multiscaling or multifractal 
properties introduced in apparently completely different systems as, for instance, models of Sclf-Organizcd-Criticality, 
DLA, random resistor networks. In this sense this approach may help to rationalize the existence of very few broad 
"universality classes" found in the scaling behaviors in very different contexts. 

This work was partially supported by the TMR Network Contract ERBFMRXCT980183 and from MURST (PRIN- 
97). 
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